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We study holographic Fermi surfaces coupled to a bosonic degree of freedom using the 
gauge/gravity correspondence. The gravity background is a charged black hole in asymptotically 
AdS spacetime. We introduce a neutral scalar field with parameters such that the system is at a 
quantum phase transition point. We further introduce a Dirac field and couple it to the scalar. 
At finite N, these fields interact in the bulk. We compute the bulk one-loop contribution to the 
boundary fermionic self-energies at various quantum phase transition points. The results would give 
an embedding of the marginal Fermi liquid model by Varma et al. into a holographic context. We 
comment on important issues that arise at low energies when one wants to take the one-loop results 
more seriously. 



I. INTRODUCTION 

A theoretical understanding of high-temperature su- 
perconductors [l| is one of the major challenges in 
physics. Cupratc superconductors exhibit many uncon- 
ventional properties. They have a generic phase diagram 
parametrized by the electron (or hole) doping x, and 
the temperature T. The superconducting phase covers a 
dome-shaped region in the phase diagram. The material 
is said to be optimally doped where the phase-transition 
temperature T c (x) is maximal. 

Above the superconducting dome (T > T c ) lies the so- 
called strange metal phase. In this phase, the material 
has unusual transport properties, e.g. the resistivity has 
a robust linear temperature dependence up to very high 
temperatures. In order to fit experimental data, Varma 
et al. proposed a phenomenological scale-invariant ansatz 
for the spin and charge fluctuations, 
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where w c is a UV cutoff (a few times smaller than the 
Fermi energy). The fermion self-energy at one- loop is 
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where x = max(|w|,T). This self-energy gives qualita- 
tively different behavior compared to the Fermi liquid 
case where ImE ~ oj 2 . Near the Fermi surface, the 
quasiparticle peak is much broader and the quasiparticle 
weight vanishes logarithmically. This system was dubbed 
the marginal Fermi liquid (MFL). 

In recent years, the gauge/gravity correspondence [2- 
H| has been applied to describe strongly coupled quan- 
tum field theories at finite charge density. The gauge 
theory is dual to a gravitational theory with negative 
cosmological constant. The gravity background has an 
"ultraviolet" boundary where it approaches anti-dc Sit- 
ter space. "Infrared" boundaries may be introduced by 
black hole horizons. Dynamics near the horizon is gener- 
ically related to low-energy phenomena in the dual gauge 
theory. Recently, holographic Fermi surfaces have been 



observed in phe nomenological ( "bottom- up" ) models 0- 
fiod (see [Til. Il2| for negative results using a "top-down" 
approach). For certain regions in the parameter space, 
these Fermi surfaces are non-Fermi liquids due to their 
strong dissipation into a locally quantum critical system 
which is described by the near-horizon AdS2 x R™ space- 
time [1, [l(| • Marginal Fermi liquids can be obtained by 
tuning the parameters of the bulk spinor field Q . A main 
question is whether marginal Fermi liquids can arise in a 
more robust way without tuning the fermion parameters. 
See [IU for a related semi-holographic work. 
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FIG. 1: Local phase diagram parametrized by the scalar mass 
m 2 and the double trace deformation k. At the intersection 
of the bifurcating and hybridized quantum critical lines lies 
the marginal critical point. 



In this paper, we study the holographic system near 
"magnetic" quantum phase transitions. The local struc- 
ture of the phase diagram is seen in FIG. [TJ The mag- 
netic _order parameter is described by a neutral scalar 
field ~ 



14l | . The boundary response function of the scalar 



takes the form of JTJ) for the MFL critical point and - 
up to logarithmic corrections - for the bifurcating criti- 
cal line. We then introduce a bulk Dirac field and couple 
it to the scalar. At finite N, the boson interacts with 
the fermion. Fermion self-energy computations reduce 
to ordinary Feynman diagram calculations with effective 
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vertices which are determined by radial integrals in the 
bulk. Assuming that the self-energy contributions from 
other fields (graviton, gauge field) are less relevant, the 
one-loop fermionic self-energy will have a marginal Fermi 
liquid form at extremely low energies. Finally, on the 
hybridized critical line we obtain marginal Fermi liquid 
and non-Fermi liquid type contributions to the bound- 
ary fermionic self-energies, depending on the spacetime 
dimensions. 

The paper is organized as follows. In Section II, we 
describe the gravity background, the phase diagram and 
the two-point functions of the fields. Section III discusses 
how the computation of the bulk loops can be reduced to 
a field theory calculation. Using these results, in Section 
IV we compute the one-loop self-energy diagrams. We 
comment on the remaining important questions in the 
discussion section. 



II. INGREDIENTS 

In this section, we briefly review the basic ingredients 
of the holographic system. Let us consider a (2+1)- 
dimcnsional rclativistic conformal field theory (CFT) 
with a gravity dual and a global U(l) symmetry. We 
turn on a finite chemical potential for the U(l), and 
put the system at finite temperature. According to the 
gauge/gravity duality, this boundary field theory can be 
described by the bulk geometry: a black hole in (3+1)- 
dimensional anti-de Sitter (AdS4) spacetime. Fcrmion 
fields in the bulk are dual to (composite) fermionic opera- 
tors on the boundary. The current in the CFT is mapped 
to a U(l) gauge field in the bulk. We further consider a 
neutral scalar field in the bulk. The expectation value of 
the dual operator will be a "magnetic" order parameter. 



zero-temperature limit, the near-horizon metric becomes 
AdS2 x R 2 . For simplicity, we rescale the coordinates, 



r -> r r, (t,x) ->• — (t,x), A -> ^A , 
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to arrive at the metric 

ds 2 =r 2 (-fdt 2 + dx 2 ) + \^-, 

r J 

The horizon is now at r = 1 and we have 
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The dimensionless temperature is given by T = 3 ^ . In 
the zero-temperature limit, the near- horizon metric be- 
comes AdS2 x R 2 where the curvature radius of AdS2 is 
i?2 = R/V&- We will use ( to denote the radial coordi- 
nate in AdS2- At small temperatures, the near- horizon 
geometry is approximately a black hole AdS2 x R 2 , 
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where Co = [2irf) 1 and T is the temperature corre- 
sponding to the r time coordinate. 



B. Matter fields 

We would like to adjust the parameters such that the 
systems is near a quantum critical point. The magnetic 
phase will be characterized by a neutral order parameter. 
The dual bulk field is a scalar with action, 



A. Black-hole geometry 

The action for a vector field Am coupled to AdS4 grav- 
ity can be written as 
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where g\ is an effective dimensionless gauge coupling. 
The equations of motion are solved by [ljj [l6| 



- -- 2 , ^ R 2 dr 2 



ds 1 = — (-fdt z + dx z ) + 

R z r z 



(4) 
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with V(0) = (0 2 + m 2 R 2 ) 
scalar will be set to zero. 

Let us now introduce a Dirac field, VP with U(l) charge 
q. The fermionic action takes the form, 
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where 77 parametrizes the coupling of the fermions to the 
magnetic order parameter. 
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This geometry is a charged black hole in AdSn with 
horizon located at r$ which is determined by f(ro) = 
0. Using the AdS/CFT prescription, [i is identified as 
the chemical potential of the boundary theory. In the 



C. Fermionic Green's function 

We are going to use the following basis for the gamma 
matrices, 
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Since in the gravity background the u> a bM spin connection 
satisfies 

^ afcM e c M r c r ab = - A T r d r In (-gg rr ) , 

we can remove the spin connection from the equations 
by introducing (^) := (-gg^) 1 / 4 ^. The bulk Dirac 
equation decouples into two equations containing the two 
two-component spinors 0i and 02 
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At the UV AdS^ boundary, any solution to the equations 
of motion for the two-component spinor can be decom- 
posed as follows, 
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The solution with ingoing boundary conditions at the 
horizon will be denoted V'in- The normalizable and non- 
normalizable solutions will be denoted ifo and 0jj , respec- 
tively. They satisfy ip^~ = and ip7 = 0. Since the 
equation of motion is real, these solutions are also real. 
The bulk Dirac equations for the two boundary spinor 
components 0x,2 have decoupled. Thus, the (diagonal- 
ized) retarded Green's function of the dual composite 
fermionic operators is given by d, [l7T - l2"oj 
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Here the label "in" indicates that the wavefunctions are 
ingoing. Note that in the boundary theory, the energy lo 
is measured from the chemical potential. 

We assume that there is a Fermi surface at a given 
momentum, i.e. at zero temperature the Green's function 
has a pole. This happens when the "source" component 
vanishes, ip^uj = 0,fcp) = 0. This equation defines the 
Fermi momentum hp. As studied in detail in the papers 
@, Ull ' the Green's function takes the following form, 
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where Qk F (w) is the AdS 2 fermionic correlation function. 
In rcscaled units, it is given by 
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Bosonic Green's function 



Generically, the low-energy behavior of the correlation 
function is dictated by the near-horizon AdS2 geometry. 
The scalar A0IS2 correlation function is Q , 
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in rcscaled coordinates. 



At low energies, the full AdS^ correlation function can 
be constructed by a matching procedure similarly to the 
fermionic case. Fixing the charge, the Green's function 
will depend on the bulk mass of the field. 

If the scalar mass is below the Brcitcnlohncr-Frecdman 
bound, then a bulk soliton interpolates between the mini- 
mum of the potential and the UV region where the scalar 
field is uncondensed. Deep in the infrared, a new AdS 2 
describes the low-energy physics (see [Til. l2lj). 

The paper [2l| introduced double trace deformations 
in this scenario. In addition to the mass, double trace 
operators can be used to drive the system into a symme- 
try breaking phase. The phase diagram is shown in FIG. 
[T] The two-dimensional parameter space is spanned by 
the scalar mass Am 2 = m? — m 2 and the coefficient k of 
the double trace deformation. 

In the following, we will briefly summarize the bosonic 
spectrum for the bifurcating and hybridized critical lines 
and for the marginal critical point. 



1. Marginal critical point 

If we set m 2 = m 2 , then solving the Klein-Gordon 
equations in the IR AdS 2 region gives the IR correlator 
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where ipo is the digamma function. At the marginal crit- 
ical point, this will be the full AdSi correlation function 
as well, see (22J. The function is analytic on the upper 
half plane, i.e. there is no condensation instability. The 
imaginary part is 
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satisfies sign(w) Im^(u;) > 0. In this approximation, the 
spectral function has no momentum dependence. Note 
the two limits, 
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Thus, in the low energy limit, the bosonic spectrum bears 
a strong resemblance to the one proposed by Varma et. 
al. [23[ for describing the anomalous transport properties 
of cuprates. This observation was also made in [2^, H3|- 
See also [23 for similar results. 
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2. Bifurcating critical point 

At a bifurcating critical point, the full UV response 
function can be built by a matching procedure, 



1 



ci - logw c + xir(w,T) 



(18) 



This form matches the numerically computed Green's 
function (see FIG. HJ. For u> ^> T, the imaginary part 
goes as Imx ~ (logw)~ . 
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where c\. 2 are real and c(fc) is a complex function of the 
momentum. In the following, we will approximate it by 
a constant c = c(k = 0). We will denote the real and 
imaginary part of this constant by d and c", respectively. 



III. FERMION SELF-ENERGY FROM BULK 
LOOP DIAGRAMS 



In this section, we compute the bulk self-energy dia- 
grams in FIG. [3] The bulk-to-bulk propagator for the 
boson (x) and for the fermion (D) are defined by 
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FIG. 2: Real and imaginary parts of the boson response func- 
tion xuv at low energies (u = -10 T. . . 10 T; T/p ~ 10 -11 ). 
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3. Hybridized critical point 

At the hybridized critical point, the T = correlation 
function takes the form plj | 
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where the differential operators act on the r\ radial co- 
ordinate. A solution can be written down [26j using nor- 
malizablc (b) and and non-normalizablc (jj) wavefunctions 
(henceforth, we are dropping insignificant constant fac- 
tors) 



D(n,r 2 ;u),k) 



®)>,a(ri;w, k) \ G a p(u, k)^\,^(r 2 ;uj, k) - ^^^^(r^ w, kj) if r x > r 2 

k)G a p(u, k) - & it0t (r 2 ;u, fc)7«^) ^b,/3(r 2 ; u, k) if r x < r 2 
I 



(20) 



Here a,/3 = 1,2 arc boundary spinor indices, 7* is a 
boundary gamma matrix and G is the boundary retarded 
Green's function for the fermion. See (2|| for a derivation 
and more details on a similar calculation. The above 
solution parallels that of a scalar 

X(n,r 2 ;uj, k) = 0|,(r>; ui, k)(p in (r < ; u, k) (21) 

where "in" indicates that the solution is infalling near 
the horizon. The ingoing solution can be written as a 
linear combination of the non-normalizablc and the nor- 
malizablc wavefunctions 



in (r;u;,fc) 



(r;u,k)+ X (u,k)<fo(r;u;,k) (22) 



where x is the boundary retarded Green's function for 
the scalar. For real w, the equations of motion are real. 



Since the j), b boundary conditions are also real, we can 
write 

lm(f> in (r;uj,k) = 4(r;w,fc)Imx(w,fc) (23) 

and thus 

Im X(ri,r 2 ;uj,k) = fain; w, k)d> b (r 2 ; w, fc)Imx(w, k) 

Note that this is an exact expression. Similarly, for the 
fermion we have [2^, [27j 

p(n,r 2 ;uj,k) = t>^ a (r 1 ;uj 1 k)p al3 (ui 1 k)ty^ l3 (r 2 ;uj,k) 

(24) 

where the bulk (p) and boundary (p a p) spectral func- 
tions are defined to be the difference between the cor- 



5 



responding retarded and advanced propagators: p = 
-i(D R - D A ) and p aP = -i(G R - G A ). 

X E 



Gt,Db D E q,D„ OA D E DbOA D E OA 

FIG. 3: Bulk loop diagrams form a geometric series when the 
fermion momentum is at k — Uf- DE(ri,r 2 ) and xb^i,^) 
are the fermion and boson bulk-to-bulk propagators, respec- 
tively. Go is the tree-level boundary fermionic Green's func- 
tion. D\, = ^\,{r) denotes the normalizable fermionic wave- 
function. 

We now want to compute the sum of diagrams in FIG. 
|3l In field theory, the diagrams form a geometric series 
and the result can be simply obtained from the one-loop 
diagram. Here we are dealing with bulk diagrams which 
do not necessarily form a geometric series. The problem 



arises because propagators between loops do not gencri- 
cally factorize: D(ri,r 2 ) ^ X{r\)Y{r 2 ). 

However, at small temperatures the bulk-to-bulk prop- 
agator (f2"U)) does factorize on the Fermi surface, 

D(n,r 2 ;uj,k) sa *i, ia (ri; u, k)G a p(u, fc)* t)/3 (r 2 ; u, k) 

(25) 

with LiJ ~ and k w kp. This occurs because G blows up 
and we can neglect the non-normalizablc piece in ((20)) . 

The factorization property of bulk-to-bulk propagators 
implies that the diagrams form a geometric series, 

G = G + GY,G + GY,GJ:G+... (26) 

where G is the resummed boundary Green's function. 
The above form of the series can be understood as follows. 
The first term is just the tree-level diagram. The second 
term is the one-loop diagram. In Euclidean signature in 
3+1 bulk dimensions, 



G aKl (ifl, fc)E reiK2 (iO, k)G K2 p{ifL, k) = 

7n~^2 dri V 9(n)dr 2 ^/ g(r 2 ) tr ^D E (n,r 2 ;iQ, - iQf m , k - q)D Eyl3 (r 2 ;in, k)xE,{r 2 ,ri)iQ! m , |g|)£> E)a ,(ri; id, k) 

where DE(r\ 1 r 2 ) and XE{ri,r 2 ) are the fermionic and bosonic Euclidean bulk-to-bulk propagators, respectively. 
flE, Q (r) and D E ^ a (r) are the bulk-to-boundary and boundary-to-bulk propagators. Finally, "tr" denotes the trace in 
bulk spinor indices, and £l' m = 2i:mT with m £ Z. At low energies and near the Fermi surface, the bulk-to-boundary 
propagators can be expressed as 

DE, a (r;ioJ, k) w $>\, il3 (r;iw,k)Gp a (iu),k) D E , a (r]iuj,k) w G a p(iu, k)^ t p(r; iu), k) (27) 

Using this we have 

■ — . C rfia , 
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{GY,G) a i3 ^T^2 ( 7TT2 dr iV '9(n)dr 2 y/ 'g(r 2 ) 



x tr m E (n, r 2 ; iO - ifl' m , k - q)^ KK2 (r 2 ;in, k)G K2 p(iQ, k)x E {r 2 , r±; iQ! m , \q\)G aK1 (id, fe)*b,«i (n; ifi, k) 

Higher loop Feynman diagrams also contain bulk-to-bulk propagators between the loops. On the Fermi surface, they 
factorize into two normalizable wavefunctions, see (|25p . and the sum of diagrams forms a geometric series. 

In order to sum the series, we compute the "half-amputated" self-energy diagram. This does not contain the 
two G Q ^(if2, k) functions from the bulk-to-boundary propagators, but does include the corresponding ^b,a{fi,2) 
normalizable spinor wavefunctions. 

S Q(3 = T Y1 J J2tt)z dri V /g ( ri " )dr2 V.9 ( r 2 ) tr (£> E (ri,r2; ifl - ift' m , k - q)^^^(r 2 ; iO, fc)x E (r 2 , n;in' m , \q\)% ia (n]ifl, k) 



in;. 

Using the spectral decomposition of the bulk-to-bulk propagators, we rewrite £ 

d ' 2 1 , /— , ,— , I p{ri,r 2 ;uj 1 ,k- q) -. Im xOi , r 2 ; u> 2 , \q\) — 



f°° dull dulo \ ^ f d 2 q / p(r\ , r 2 \ wi , k — q) , 

V aB (in,k = / _A_»t£ / « driVadrav^tr p \ 2 ,' /' ^-% >p (r 2 ;iU,k) 



* bia (riiin,fc) 



From eqn. ((24]), the bulk-to-bulk spectral function p can be expressed using the one in the boundary theory, 

f 00 du>i du) 2 x — . r d 2 q , , , 

^(^ fe )=y_ oo — — te y ^ d n^ 2 v^x ( 2 8) 

. ^ q)P<T lt r 2 (ui u k~qj^ b (r 2 ;ui,k-g) 7 <j> b {n;u 2 , \q\)<p b (r 2 ;u 2 , |g|)Im X (^2, \q\) ^ . ^ 

tr — — ^t, p(r 2 ;in,k) — * b a (n; ill, k) 

\ lil — III' — Ull III' — U>2 I 



The integrand has factorized into n and r 2 dependent parts. Evaluating the sum over fl' m gives 



v m ft [°° dLo 1( ko 2 [ d 2 q f ^ A tanh(§,) +coth(^) 
7_oo 2?r 2tt 7 (27r)^ 7 S2 - uji - uj 2 + le 

Xlmx(w2, kl) tr % i<Tl (n;a;i,fc - q)p ai a 2 (uJi,k - q)%,a 2 (r2',^i,k - q)^ b ^(r 2 ;ifl, ^^(rx; ifi, £;) (29) 

Let us now pick a basis where p a p(uj, k) is diagonal and fix the two basis vectors vE f b,a=i,2( r ; k) for the fc direction. 
For other momentum vectors in the x — y plane we define a basis for the normalizablc wavefunctions 

(30) 



where 1Z and R are the boundary vector and spinor rotations, respectively. Note that we have included a phase in 
order to cancel the monodromy R(2n) = —Id. 

Since k — q is not parallel to fc, the spectral density p a p(k — q) in (|29[) has to be diagonalizcd, 



*b,a(»"i;wi,£- q)p a p(ui,k- q)^^p(r 2l u}i,k- q) = (e * i? QCTl >J b CTi J R a<J2 p a2Cr3 (k - q) (R asP e^R^Jk 



p a p(k\k-q\) 



where k = k/\k\, and ip is the angle between k — q and k. The rotated wavefunctions can be expressed using the 
original ones in the k direction 



e £R a p% t p (ri\u,k - q ) = ^ b . a [ n; u, k\k - q\ 

Then, we can rewrite (|29|) as 

diO\ du> 2 f d 2 q 



(31) 



T, a p(iQ, k) = 8 af3 



tanh (§£ +coth H 
dr 1 ^gldr 2 ^gi — ft^^ — , <ftb(n;^2, |s|)0b(^; w 2 , |g|) x 



27r 27r 7 (27r) 2 J -"-v^j- i v ^ fl — uji — u! 2 + ie 
xlmx(w 2 , \q\) tr *b,a(n; wi, fe|A: - g|)p QQ (wi, - gl)*^^; wi, fe|A; - g|)*b,a(r 2 ; iO, fc)* bja (ri; iQ, fc) (32) 
There is no summing over a. The momentum arguments inside the trace are now all parallel vectors. 



A. The effective vertex 



We will be concerned with the boundary component a for which there is a Fermi surface, i.e. G a ^(0,kp ) = at 
T = 0. The corresponding self-energy in (|32|) can be rewritten (we will suppress k in the following) 

E aa (Q,k) = ^2 [ "^"^|2^p" Ab ' Q w i'|£ _< z1; w 2,|g|) A bi « (wi,|fc-g|; Q, \k\; UJ 2l \q\] x 



xpoofwi, |fc - g|)Imx(w2, M) 



tanh(^) +coth(§|) 



(33) 



Q — u>i — uj 2 + ie 

where we defined the effective vertex using normalizablc two-component wavefunctions 

Ab, Q (wi,£i; w 2 ,fc 2 ; w 3 ,fc 3 ) = y rfr % /g7/i b Q (r;wi,A:i)V'b, a (f;w 2 ,A: 2 )(?!)b(r;w3,4) 

We will use the 2x2 gamma matrices 7- = ztr 1 , 7- = u 3 , 7— = — a 2 . In order to compute Im S, we do the tu 2 integral 

du x d d ~ x q 



Im £(fi,jfe) 



2tt (2tt) 

X/3 aa (wi, |fc - g|)Imx(fi - wi, |g|) 



2 Ab, Q (0, |fc|; wi,|fc-g|; fi-wi,|gl) A b Q (wi,|fc-gl; n,|ft|; 0-wi,|gl) x 



tanh \^^;J + coth 



2T 



(34) 
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At low temperatures, the hyperbolic functions sup- 
press the integrand outside the window ui\ = [0 ... £1} . 
The fermion spectral function can be approximated by 
ImGfl w Z5(k± — Re£ ). This effectively restricts the 
momentum integral to the vicinity of the Fermi surface 
and we have |fc| ~ k F , \k — q\ w k F . 



1. Infrared divergences: T £ 1 <C a; <C /i 

A potential divergence may come from the deep IR 
region: T <C C 1 "C w. (See Q for divergences in a 
similar effective vertex in AdS n ). Let us fix T = and 

use the AdS 2 metric: ds 2 = ~ dT ^f dc - . The neutral bo- 
son and charged fermion wavefunctions can be expressed 
using Hankel and Whittaker functions, respectively [29| 



^(C)=c out H«K) + c in H(, 2 )K) 



(35) 
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<KC) = ^=(1 + ia 1 )^ = c out W_ 4 _ igy (2ia;C) 



fh + im 



where the notation a 3 in the index of the Whittaker func- 
tion indicates ±1 when acting on the top/bottom com- 
ponent of the spinor, and in our units fh = — ( — l) a k. 

Let us first look at the effective vertex with ingo- 
ing/outgoing boundary conditions for the wavefunctions. 
For uj^ 3> 1, the ingoing wavefunctions 

</£(C; w) ~ c^e^C 0in(C; w) ~ CV" C 

where the two signs correspond to the two spinor compo- 
nents. The outgoing wavefunctions satisfy i/'out = (V'in) j 

0out = (</>in)*, 

fe(C;w) ~ Fle-^C* 0out(C;w) ~ C°e-^ c 

The Yukawa coupling can be written as 



V2 



V2 



From the last form we see that the Yukawa coupling com- 
bines ip + with ip~ . For ingoing wavefunctions this gives 

V'inV'in^in ~ i ("0it) <fe<fcn + i (^in) V'in^in ~ C° 

Note that the ( ±lq factors have canceled each other. 
There is no divergent contribution to the effective ver- 
tex 



cutoff S 



A in (Q,0Ji,UJ2,k F ) = j TJ-V'inV'in^m < OO 



Similar results apply to outgoing wavefunctions. Let us 
now turn to normalizablc wavefunctions 



Ipin - 4>ou 

2ilm G 



(36) 



The Yukawa coupling contains various terms 

The last two terms couple the spinor components ip^ to 
■0~ ut , which scales for large ( as 



(37) 



As C — > oo , this is more divergent than the purely ingoing 
case. The effective vertex has a logarithmic divergence 
at q = 0. For positive fermion charges, 



d( 
T 



,i(0+ie)C fliq 



C 21 " < OO 



Note that the finiteness of the vertex crucially depends 
on the spinor index structure and is not guaranteed for 
other types of interactions. 

The UV region (i.e. the bulk of AdS 4 - RN) con- 
tributes a constant to the effective vertex. We have seen 
that the contribution from the region C" 1 <C ui is not di- 
vergent. It remains to study potential divergences from 
the region ui <C C _1 - 



2. Infrared divergences: T, uj <C C 1 ^ A 1 

In the bulk, ip\, and <p\, are temperature-independent 
to first order. These functions match onto those in the 
near-horizon AdS 2 region: 



( | 



( |f 



where C = —^-r is the rescalcd near-horizon coordinate. 

*5 r— 1 

The IR contribution then scales as 

phorizon 

A Jfl (0,0,fc F ) = / drv^^(r)Vb(r)0b(r) 

J cutoff 



rp / £ \ 1—2uf / 



rp2v F —^+V B 



Here T c is a constant which sets the cutoff that separates 
the IR and UV regions. Thus, the contribution from the 
IR AdS 2 is not singular in the zero temperature limit if 

Note that when the internal fermion propagator is not 
on the Fermi surface, the corresponding wavefunction 

V>jf (C;0, fc) has a piece that goes like ( . The ef- 

fective vertex then scales as T VB ~^ which is divergent 
when vb <\- However, the full expression for ImE con- 
tains a factor Im Gr ~ x 2l,p which - for large enough v F 
- renders this contribution harmless. 
Finally, we can substitute 

A(fi, wi, Q - wi, k F ) -> A(0, 0, 0, k F ) = A(k F ) (38) 
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The resulting expression is now familiar from finite tem- 
perature quantum field theories 



B. Bifurcating critical point 



Im S(fi,fc) « K{k F f 



^ll_|lmG R (a;i,|fc-gl) function is [2 
2ir (27r)^ 



ximx(n-wi,kl) 



tanh 



^1 \ 
27V 



coth 



£7 — Wj 
2T 



(39) 



Let us consider the T«w limit. The bosonic spectral 



Imx(w) 



l0 S £ + h 



(46) 



The only consequence of the radial direction is that we 
had to replace the interaction by the effective vertex A. 



IV. FERMION SELF-ENERGY RESULTS 

In this section we will use the previous results to com- 
pute the one-loop bosonic contribution to the fermion 
self-energy in various regions of the phase diagram FIG. 
[TJ As it is clear from eqn. (f3~4"|). this will be an ordinary 
field theory calculation where the interaction vertices are 
determined by radial integrals in the bulk. 



This is again independent of the momentum. At T = 0, 



the factor tanh ^ +coth 
wi e (0...Q). Thus, 

ImS(n) r 



2T 



has support in the interval 



duii 



1 



2tt log 2 

At large ui a and il <C ui a we get, 
ImS(n) 



n 



lo: 



,2 n 



(47) 



(48) 



A. Marginal critical point 



C. Hybridized critical point in 2d 



The boson spectral function assumes the momentum- 
independent scale-invariant form, 



Im k, T) = tanh 



(40) 



This approximation is valid for small momenta. Near the 
Fermi surface the fermionic spectral function becomes a 
Dirac delta, 

lmG R . (k F + k ± ,ui) « Z5{k ± - Re E (w)) (41) 

It is easy to evaluate the momentum integral, since only 
Gb.,o contains q. We get, 



ImS 



en) ~ r ** 



tanh tanh h coth 

2T V 2T 2r 



After performing the ui\ integral we get 



ImSffi) - ft coth — 
2T 



(42) 



(43) 



The real part can be computed using the Kramers- 
Kronig relation. The full result is, 



S(w,T) 



inT 
u 



log 



2^f 



2ttT 



(44) 



where denotes the digamma function. The term in- 
volving the UV scale ui c renders the large frequency limit 
well-defined, i.e. temperature independent. This func- 
tion is a smooth version of the original marginal Fermi 
liquid self-energy: 



£mfl(w,T) ~ wlog 



7T 

i— x 
2 



(45) 



where x = max(w, T). 



In this case, we have to take into account the momen- 
tum dependence of the susceptibility. The fermionic self- 
energy, 



£(fi,fc) « K{k F f 



2 f°° duii du)2 f d 2 q 



2ir 2ir J (2tt) 2 



(49) 



xImG R , (wi, \k - q\)lmx(ui 2 , \q\) 



tanh(i%) +coth(^; 



Q — u>\ — ui2 + ie 
At T = 0, the imaginary part becomes 



Im£(Q,fc) - sgnftA(fc F ) 



dwi f d 2 q 



lo 2tt 7 (2tt) 2 
xlmG Ri0 («i, \k- g|)Imx(f2 - w ls |g|) 




(50) 



FIG. 4: Variables in momentum space. 

The Dirac delta in Im G restricts the momentum inte- 
gration such that \k — q\ — kp = Let us ignore the 
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curvature of the Fermi surface and denote the tangential 
component of the momentum by x (see FIG. HJ) - Then, 



ImS(fc_L,r2) - sgn^A(7c F ) 2 

x J "2tT J lm x(^,q) 



(51) 



2. v B = l 

Let us neglect logarithms and approximate the suscep- 
tibility by the simple function, 



oj— Q— loi; q—\ I 



We will need the imaginary part of the susceptibility, 



Im^(u, k) 



(Clfc 2 + C2W 2 + c'w 2 ^) 2 + C" 2 UJ^B 

In the following, we will assume 51 > 0. 



(52) 



E" ~ ai 1_r \ 


I E" ~ oj 2t " 2 


r ~w 2v+l \k^\- 3 



0.0 0.5 



1.0 

v 



1.5 2.0 







FIG. 5: Scaling of the one-loop fermion self-energy 
ImS(u, fcx) along the hybridized quantum critical line in 
(2+l)-d. v = vb denotes the bosonic scaling constant. 



1. VB > 1 

The imaginary part is approximated by 
Im^fw, k) 



(53) 



(ciA; 2 + c 2 uj 2 ) 
We obtain the following limits, 

for < fi < \v F kj_ 



ImE~ 



2c{(l+2u B ) ' IfclP - 



^j— for Q = \vpkj_\ 



{ 4j^-^- 2 for \v F k ± \«n 



(54) 



/ 3 3 / co /Co 

/ = ir x l + 2i/ s , -,-,2 + 2v B ; 1 + —v F , 1 - iJ — v F 
\ 2 2 V c i V c i 

where f\ is the Appell hypergeometric function of two 
variables. There is a crossover at the energy scale set by 
vpk±. Above this scale, the self-energy becomes inde- 
pendent of the momentum. 



x(w,fc) 



1 



C\k 2 + CLO 2 



(55) 



The fermion self-energy is then a scaling function 
Im E(fi, k) = f ■ In the two limits, 

f ^rfir for \n\«\v F k^\ 
ImE~ I (56) 

[c for |n|»|t>Ffc±| 

where the constant is 



C : 



c"nai l)VTTa + alog [(\A + £- X ) (vT+a-1 



with a = 

Cl 



2c%v F (l + a): 



3. < I/ S < 1 



The imaginary part of the susceptibility can be ap- 
proximated by 



Imx(w,fc) 



(Clfc 2 +C'W 2 ^)^ +C" 2 LJ^B 



(57) 



The transverse momentum integration is straightforward, 
and the frequency integral can be done in two limits 1 ' 2 



° 2 " s+1 for |n| <£a\kx\ 1/vB 



ImS 



where a 



2c((l+2u B ) ' \k±\ d 

for a\k A l 1 /^ « Ifi 



/ci|c| 



1 — VB 



\ci/c\- 



and 



c = ce 



(58) 
Near the 



marginal critical point (vb ~^ 0), the small fre- 
quency region rapidly disappears and the self-energy be- 
comes momentum- independent. Note also that the low- 
frequency regime has the same form as that of the vb > 1 
case. 



D. Hybridized critical point in 3d 

Here we assume that the results of Section III apply to 
3+1 boundary spacctimc dimensions. Then, similarly to 



1 This form of the self-energy has been confirmed by numerical 
calculations and an analytical calculation at ug = i where the 
integrals can be evaluated exactly. 

2 We found the same result in |30H ; see eqn. 136 on page 1040. 
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the (2+l)-dimcnsional case, we can write 

f n d 

ImE(fci,l]) ~ sgnftA(fc F ) 2 / — 

Jo 2lT 



(59) 



/ dx 27nrlmx(w,q) 
Jo 



w=fi— wi; ^r)" 



Compared to the 2d case, the integrand has a 27rx factor 
which comes from the momentum integral in the extra 
dimension. 




0.0 0.5 



FIG. 6: Scaling of the one-loop fermion self-energy 
ImE(cd,fcx) along the hybridized quantum critical line in 
(3+l)-d. v = vb denotes the bosonic scaling constant. 



1. v B > 1 



We have, 



c"tt n J "-sH 



c{(l+2u B ) |/=_l| 2 



for < 9, < |w f ^_l| 



Im £ - ^ 



■ • for = \v F k±\ 



s^l.^B-i for j^fcj < ft < f2 uv 

C l C 2 

(60) 

where / can be expressed in terms of a hypergcomctric 
function of vb , Ci 2 an d w f ■ 



Im £ 




° 2 " B+1 for |J2| < al/feil 1 /"* 



for Qflfcxl 1 /^ < |Q| < Q uv 
(62) 

E. Away from the critical point 

In the bifurcating case, the boson spectral function has 
the frequency dependence, 



XM ~ cuj 



2v B 



(63) 



where c € C. We neglect the momentum dependence of 
Vb- Then, 

ImS(O) ~ r^ll mc (0- Wl f B ~Im (cft 2 ^ +1 ) 

When v B < 0, this locks the frequency scaling to be lin- 
ear. Although the formula suggests a log-oscillating self- 
energy, it can only be trusted at high enough frequencies 
where the first oscillation has not happened yet. Below 
this energy scale, the scalar soliton kills the oscillating 
behavior. 

Near the hybridized critical line, 



A;) 



1 



S + cik 2 + c 2 uj 2 + c(k)ui 2vB 



(64) 



where S parametrizes the distance to the QCP. The result 
for vb > 1 is 



Im£(ft,fc) 



Im (cfl 2vB+1 " 



6 + max{£,k±y 



DISCUSSION 



2. vb — 1 

The self- energy, 

f §^ for M<\v F k ± \ 
Im£~| (61) 
[ S-Q, for |0| > \v F k±\ 

where the constant is 



C = 7r / dw arctan 

2 Jo 



In this paper, we studied the holographic response 
function of fermion operators at various quantum phase 
transitions. The order parameter is a neutral scalar field 
in the bulk which may represent antiferromagnetic order. 
The phase diagram is shown in FIG. [TJ 

When the scalar m 2 reaches the AdS 2 Breitenlohner- 
Freedman bound at the 'bifurcating' critical line, a clas- 
sical instability develops and the system enters the or- 
dered phase. Certain bulk loop corrections to the bound- 
ary fermionic two-point functions can be reduced to or- 
dinary field theory calculations where the vertices are 
replaced by effective interaction vertices. These effec- 
tive vertices are determined by radial integrals in the 
bulk. At the 'marginal' critical point, the one-loop bulk 
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contribution from the scalar to the fcrmionic self-energy 
is £(!) ~ a; logo;. If the absence of other relevant low- 
energy effects, the Green's function can be written as 



G(w,k) 



hi 



(fc - k F ) 



J (o) 



N- 2 E (1 



(i) 



(65) 



which would give a marginal Fermi liquid. Note that 
bulk loop diagrams arc suppressed by powers of N and 
the one-loop results can only be important at very small 
energy scales. In the case of a marginal Fermi liquid, this 
may happen at exponentially low energies where the log- 
arithm cancels the N~ 2 factor. However, at such energy 
scales it is important to take into account the backreac- 
tion of the bulk charge density on the background geom- 
etry. This deforms the near-horizon AdS2 into a Lifshitz 
space [3lM35] and may invalidate the one-loop result. 

Double trace operators can be used to drive the sys- 
tem into the symmetry breaking phase through the 'hy- 
bridized' quantum critical point. The one-loop results 
in two spatial boundary dimensions give En) ~ lu 1 ^ Ub . 
Here fe(m 2 ) is a bosonic scaling exponent which van- 
ishes at the marginal critical point. In this case, it is im- 
portant to take into account the Landau damping of the 



bosons. This presumably changes the low-energy physics 
in important ways. 

The effective Yukawa vertex in this paper was finite. 
The finiteness crucially depends on the spinor index 
structure. It would be interesting to study other bulk 
interactions where the effective vertex has an infrared 
divergence (see also [H])- 

Finally, it would be interesting to study a similar sys- 
tem using a consistent 'top-down' approach where the 
boundary field theory and the operators corresponding 
to the bulk fields are explicitly known. 
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